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Abstract 

Fluid transport in microfluidic systems typically is laminar due to the low Reynolds number 
characteristic of the flow. The inclusion of suspended polymers imparts elasticity to fluids, allowing 
instabilities to be excited when substantial polymer stretching occurs. For high molecular weight 
polymer chains we find that flow velocities achievable by standard electro-osmotic pumping are 
sufficient to excite extensional instabilities in dilute polymer solutions. We observe a dependence 
in measured fluctuations on polymer concentration which plateaus at a threshold corresponding 
to the onset of significant molecular crowding in macromolecular solutions; plateauing occurs well 
below the overlap concentration. Our results show that electro-osmotic flows of complex fluids are 
disturbed from the steady regime, suggesting potential for enhanced mixing and requiring care in 
modeling the flow of complex liquids such as biopolymer suspensions. 
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I. INTRODUCTION 



It is known that the inclusion of polymers into fluids allows flow fields and instabili- 
ties associated with elastic effects to be excited in pressure or other mechanically-driven 
liquids pp. Here we show that by adding dilute amounts of inert polymer to liquids exten- 
sional instabilities can be excited in microchannels under steady electro-osmotic flow (EOF) 
pumping. 

For mechanically-driven liquids there is always a significant shear component to the fluid 
flow, due to no-slip boundary conditions; the effective slip boundary conditions of EOF sup- 
presses the shear allowing elongational effects to be emphasized and studied. Considerable 
efforts have considered "infinitely dilute" single-chain behavior of DNA in elongation flow 
(see, for example, Refs. [2J, [3] ) , making use of electro-phoresis which places a body force 
on the polymer coil dragging it through a stagnant liquid. We consider here finitely dilute 
solutions driven at the boundary and undergoing elongtional stress in order to investigate 
the possibility of instabilities in such flows, which we determine can be easily excited for low 
viscosity solvents. 

Placing polar liquids in contact with ionizable surfaces produces an electric double layer 
(EDL) where ions of opposite charge are attracted to the wall [3]. This electrically non- 
neutral EDL sheaths surfaces, allowing electric fields to produce flows by driving the sheath, 
which in turn viscously couples to the bulk fluid away from the wall [3]. Such electro- 
osmotic flow is desirable for use in micro- and nano-fluidics, as flow with flat velocity fronts 
is observed away from the EDL in uniform channels of arbitrary cross section, with shear 
confined to the EDL regions. As a result scalability is superior to pressure-driven flows 
which become difficult to drive through very narrow channels. 

Microfluidic devices usually display creeping laminar flows [5], for both pressure and 
electro-osmotic pumping, as the achievable Reynolds numbers are small when working with 
small volumes in microscaled (0(100) //m typical) channels. The low Reynolds regime 
prevents instabilities and turbulence which enhance mixing of chemical reactants allowing 
timely and complete reactions. As counter examples to typical "smooth" low Reynolds 
number flows see, for example, Ref. [6] where high velocities result in intermediate Reynolds 
numbers (0(1-100)) allowing inertial effects to induce Lagrangian chaos and efficient mixing 
in planar 2D microchannels and Ref. [7j where bas-relief (3D microchannels, requiring more 
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demanding microfabrication) of a staggered herringbone shape allows Lagrangian chaos and 
effective mixing for Re < f . 

Polymer coils in solution can be stretched in flows. Thermal agitation causes a relaxation 
back towards the globule coil state and polymers act as entropic springs [lj. Polymeric 
solutions therefore have a relaxation time, A, and if flow times, £f = Lchar/v, where Lchar 
is a characteristic channel dimension and v is the flow velocity, are driven to be shorter 
than the relaxation time A significant elastic stress can develop and lead to secondary flows 
and instabilities (see Chapter 2 of Ref. pQ for a general overview). The Deborah number, 
De — A/if, is a dimensionless parameter describing viscoelastic solutions; when this number 
reaches a geometry-dependent threshold, roughly of order unity, instabilities often occur [1] . 
For example, it has been shown that chaotic instabilities can arise in Dean flow of polymer 
solutions through curvilinear serpentine microchannels [8J, where earlier characterization of 
macroscopic-equivalent flows [9] demonstrated the presence of elastic turbulence [10J. In 
addition to these flows were elastic turbulence and chaotic flow are established it is known 
that rapid constriction/expansions can also lead to instabilities pQ, for example modifica- 
tion of serpentine microchannels |8J to add asymmetric corrugation (triangular constrictions 
followed by sudden expansions along the channel) demonstrates the formation of unstable 
vortices [11] when polymeric fluid is pumped through the channel at sufficiently high rates. 
These unstable vortices are excited in linear channels as well and display apparent random 
growth and decay dynamics [12] which, due to the corrugation asymmetry, display differing 
dynamics depending on flow direction. The unstable vortices arise at Re < 1 for highly 
elastic polymeric liquids; in the same device Newtonian flow (water) also displays vortices 
at intermediate Reynolds numbers (Re ~ 5) due to inertial effects. 

II. EXPERIMENTAL SETUP 

Microchannels were formed in glass with 100 /mi and 200 /mi width and 20 /an depth 
to the bottom of the "D"-shaped chemically etched cross sections. Forming 2:1 constric- 
tions (see Fig. [l]a and b) causes flow acceleration and deceleration as liquid transverses 
through the constricted region, inducing an extensional (stretching) flow field. Polymer so- 
lutions were made by adding 18 x 10 6 Da high molecular weight (HMW) or 5 x 10 6 Da low 
molecular weight (LMW) polyacrylamide (PAAm) from Polysciences Inc. to a 20:80 vol.% 
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methanol:water mixture, where methanol is used to prevent aging (chain scission) during 
storage |T3]. Neighboring polymer coils impinge on each other at the overlap concentration 
c* of ~ 300 ppm and ~ 850 ppm for the HMW and LMW samples, respectively. We find 
concentrations both below and mildly above the overlap concentration of the HMW sample 
display similar behavior (as seen in Fig. [5]b; see also Fig. [6]). 

The glass microchips are loaded into a Microfluidic Toolkit [50] and flow is electro- 
osmotically driven with voltages of up to 0.4 kV (electric fields as high as ~ 900 V/cm). 
The interface between the two input streams is observed by fluorescently tagging one input 
solution with a small amount (« 21 fiM) of tetramethylrhodamine (TAMRA) dye and ex- 
citing with a 532 nm (green) diode laser. Laser-induced fluorescence profiles are collected 
for various flow rates by focusing the laser to l/10th of the channel width (10 /im beam 
waist) close to the dyed-undyed fluid interface centered in an extensional region (focal spot 
indicated with a small circle in Fig. [Tjb) , and collecting light emitted from the dye at 570 
nm (yellow) with a photomultiplier tube (PMT) after filtering with a 568 nm bandpass filter 
(10 nm FWHM). 

Data runs were 100 s in duration, with the first 20 s removed prior to analysis to allow 
startup and transient behavior to attenuate. Care was taken to thoroughly clean and rinse 
microchannels before use to ensure clean surfaces to prevent pH shifts or other effects from 
modifying flows. The cleaning protocol was to flush the channels with: water (x2), ~ 1 
M nitric acid, w 1 M NaOH, water (x2); first by filling all the reservoirs and then pulling 
solutions through the channels by applying vacuum to the output reservoir (reservoir 4 in 
Fig. 0a). 

Flow velocities in EOF are given by v = ^eofE, where /aeof is the electro-osmotic 
mobility and E the applied electric field. By timing the movement of a marker of dye over a 
known distance the electro-osmotic mobility was measured to be (5.6 ± 0.5) x 10~ 4 cm 2 /Vs 
for polymer-free solutions, which agrees with the value of 6.00 x 10 -4 cm 2 /Vs reported in 
50 /im diameter fused silica capillaries [IB] . The HMW solution at 120 ppm has a /xeof of 
(5.7 ±0.7) x 10~ 4 cm 2 /Vs, indicating that for dilute solutions of large polymers the electro- 
osmotic mobility is not significantly modified from the pure solvent case. For the LMW 
solution at 120 ppm the mobility is reduced to (3.1 ± 0.2) x 10 -4 cm 2 /Vs, possibly due to 
polymer coils interacting with the Debye layer [T7] which is on the order of 100 nm thick 
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FIG. 1: Device schematic. The microfluidic device (a) consists of four-reservoirs (two input 
reservoirs at positive voltage that flow into the grounded fourth reservoir, the third reservoir floats 
during experiments described here) connected via micro channels. A dyed (reservoir 1) stream 
merges with an undyed (reservoir 2) stream and flows towards the ground (reservoir 4) through 
microchannels defined by a 20 ^m deep chemical etch which results in the distinctive D-scooped 
profile associated with wet etching. Past the intersection that merges the two input streams 200:100 
fim wide channel constrictions occur creating extensional flow fields. The small (not to scale) circle 
in (b) indicates the photomultiplier measurement spot (~ 10 /im focused waist of a 532 nm diode 
laser) used to measure fluorescence emitted by excited dye, while the large dashed circle outlines 
the field of view for microphotographs. 

We found that use of a LMW sample from Polysciences Inc. which contained ~ 1% 
ammonium chloride stopped electro-osmotic flow, even after dialysis. As PAAm is known 
not to adhere to channels [19J, and ~ 1% ammonium chloride solutions in water:methanol 
solvent demonstrated normal flow behavior, the mechanism preventing flow in the sample 
which included salt is not known. One possibility is the presence of a contaminant copolymer 
similar to poly-N-hydroxyethylacrylamide [20J, which adheres to side- walls and significantly 
reduces the electro-osmotic mobility. 
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III. INSTABILITIES IN EXTENSIONAL MICROFLOWS OF POLYMERIC SO- 
LUTIONS 



Creating a visually observable fluid interface by flowing dyed fluid from reservoir 1 and 
undyed fluid from reservoir 2 to reservoir 4, Fig. [TJa, allows us to investigate fluid flow 
properties. The channel geometry is designed with 2:1 constrictions in order to induce 
extensional flow and the fluid interface is monitored in the center of this region (at the 
white spot in Fig. |TJb) . The Reynolds numbers here are < 0.01, below the nominal creeping 
flow limit of Re < 1. Without polymer creeping laminar flow persists to the highest applied 
voltages (Fig. |2}b), with the PMT traces appearing similar to the no-flow case (see Fig. |2}a). 

(a) 
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FIG. 2: Normalized photomultiplier tube time traces. A PMT monitors the interface of the 
dyed/undyed streams of methanol- water solutions with (a) no driving voltage (no flow), (b) driven 
at 0.4 kV, and (c) driven at 0.2 kV with 120 ppm of HMW polymer added. With no polymer the 
PMT trace of the flow is smooth (b), albeit with noise from inherent laser diode instability as also 
seen for no flow (a), indicative of a stable interface and the laminar flow typical in microchannels. 
With polymer added rapid interface motion and instabilities are observed at moderate driving 
voltages, as evident by the dramatic fluctuations in the PMT signal (c). 



In contrast, very large fluctuations arise, at sufficiently large driving electric fields, when 
polymer is added to the solution (Fig. [2jc). This arises from an instability of the flow causing 
the local dye concentration to vary under the green excitation spot. The presence of polymer 
results in large movement of the dye/undyed interface relative to the illumination spot, 
generating order-unity fluctuations in Fig. [2jc. We note that extensional flow also exists in 
the cross-channel region, due to acceleration, which excites instabilities in polymeric fluids 
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at high flow rates; the Deborah number in the cross-channel region is lower than in the 
designed constriction and these instabilities are eliminated at low flow rates. 
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FIG. 3: Autocorrelations of photomultiplier tube signals (Fig. [2] data). Methanol- water solutions 
with (a) no driving voltage (no flow), (b) driven at 0.4 kV, and (c) driven at 0.2 kV with 120 
ppm of HMW polymer added. The autocorrelations of polymer free fluid flow decorrelate within a 
few sampling time steps and no substructure is apparent, while unstable polymer flows show rapid 
drop off followed by substructure without prominent peaks, indicating the presence of a multitude 
of frequencies. 

Photomicrographs (Fig. [4]) captured by wide field excitation of the dye using an external 
laser show the stable interface found in polymer-free flows (Fig. |4ja), sampling this interface 
with a focused laser spot and a PMT results in a smooth time trace (see Fig. |2]b). Adding 
polymer leads to fluctuations, of which the largest scale features are apparent in Fig. |4jb and 
c (compare with Fig. |2}c). The largest spatial scale of viscoelastic instabilities is determined 
by the system size, as for inertial instabilities pp. There is no clear lower bound on spatial 
scale, with polymer size, separation, or diffusion lengths being possible cutoff lengths |2T] . 

No dominant oscillations were observed in PMT signals collected. The autocorrelation of 
PMT traces of polymer solution flows have quick decay times (~ 0.1 s decorrelation time, 
taken as the e-folding time), and displayed substructure indicative of multiple frequency 
components while not exhibiting sharp, large, repeated ringing which denote dominant si- 
nusoidal frequencies. Autocorrelations of PMT signals for polymer-free fluid decorrelated 
in « 2 sampling time steps and showed no substructure. Fig. [3] shows the autocorrelation 
function for Fig. [2] data. Scatter and spectral plots (not shown) also did not show evidence 
of distinct sinusoidal frequencies [23] in unstable polymer flow. Taken together these mea- 
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sures indicate that the flow is chaotic [8]. More detailed evidence is required in order to 
determine if the system investigated here displays true elastic turbulence [21], for example 
measurement of flow velocity fluctuations and their spectral scaling properties. 
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FIG. 4: Photomicrographs of (a) no polymer and (b, c) 120 ppm HMW polymer flows at 0.2 
kV. Image (c) is taken one second after (b). Dashed lines indicate field of view and microchannel 
constriction outline. The interface between dyed and undyed liquid streams is (a) stable for polymer 
free samples and becomes (b, c) unstable with the inclusion of polymer. The scale bar is 100 /im. 

The PMT signal is dependent on dye concentration, possible thinning effects [21] , and flow 
velocities. To reduce such effects the PMT signal is normalized by the average voltage. To 
study the development of the instability versus driving voltage, the PMT signal fluctuations, 
5pmt, defined as var f(^w> where var(V) is the variance of the PMT voltage signal from 
the average (V), are measured for various driving potentials. The precise location of the 
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excitation spot from the interface will affect our measure 5pmt- Our absolute positional 



accuracy is roughly 5 /mi and, by measuring away from the interface on either side for 
a given experimental alignment, we estimate the introduced error to be roughly 10% of 
the measured fluctuation for unstable flows at higher flow rates (0.2 kV and above driving 
voltages) which is also the estimated error found by comparing differing experimental runs 
(which requires realignment). At lower flow velocities the system noise begins to dominate, 
and for no- flow conditions run-to- run differences in measured fluctuations is ca. 30%. 

Polymer-free solutions display laminar flow for all driving voltages (Fig. [5ja). In contrast 
to this laminar flow, fluctuations grow in size as the driving voltage increases for a concen- 
tration of 120 ppm of the HMW sample (Fig. [5ja). As seen in Fig. [5}b similar behavior 
occurrs at other concentrations of the HMW sample, both below (60, 90 ppm), and above 
(480 ppm) the overlap concentration. 

Scaling properties of viscoelastic instabilities are not well understood, but the basic ele- 
ments of curvature and elastic normal stress (stretching) along streamlines are well estab- 
lished as underlying instability [25] as these mechanisms lead to "hoop-stress" pQ allowing 
stretched polymers to cross streamlines, disturbing laminar flow. Fluctuations should grow 
with concentration (as the movement of individual polymer coils underlies the fluctuations). 
Changing the concentration and investigating fluctuations at a constant (0.2 kV) voltage 
indicates that fluctuations rapidly grow from the noise floor and then peak and plateau with 
a gentle decrease as we increase concentration, with fairly uniform measured fluctuations 
from p:X to ~ IX the overlap concentration, Fig. p| 



Our fluctuation measure var takes its maximum possible value of 1 for square- wave 

signals; the precise maximal value for a given signal is < 1 and depends on the nature of the 
waveform. Smooth growth in our measure occurs until the maximum, for a given waveform 
with linear amplitude growth. When the maximal value is reached a sudden apparent 
threshold will occur and our measure will saturate. The maximal threshold observed here 
is followed by a mild decrease rather than a strict capping and saturation; this is suggestive 
of a qualitative change in behavior rather than an artifact of our measure. 

Modeling the polymer solution as a colloidal solution |26j one expects hydrodynamic in- 
teraction mediated molecular crowding effects above a volume fraction cf> of approximately 
0.1 [27] where such effects are empirically found [51J. To estimate the concentration when 
hydrodynamic interactions between polymers play a significant role we take Cthresh = hh^ssh^ 
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with (pthresh empirically known [27] to be ~ 0.1 and 0* corresponding to the overlap con- 
centration which, depending on the particular packing taken near overlap, ranges from 
roughly 0.5-0.7; this leads to an estimate of c t hresh ~ 40-60 ppm for our HMW polymer; 
this estimated threshold is close in value to the observe peak in Fig. [6j The similarity in 
thresholds suggests that crowding may be the mechanism inhibiting growth in measures 
fluctuations; the tendency for polymer coils to deform makes the hard colloid model, which 
assumes polymer coils act as hard spheres, only approximate. It is known that molecular 
crowding strongly affects and reduces molecular transport [2E] by hindering movement of 
macromoleclues; hindered movement of polymers across streamlines will limit fluctuations, 
as we observe here. 
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FIG. 5: Fluctuation amplitude in normalized fluorescence signals. Fluctuations in the motion 
of the dyed/undyed fluid interface are measured with normalized PMT voltage for (a) 120 ppm 
of HMW polymer; the error is roughly on the order of 10% at higher flow rates for unstable 
flows. Control (polymer-free) solution fluctuations are nominally zero (compare with Fig. [2]b). 
The instability threshold is predicted to be ~ 0.02 kV. Similar behavior is displayed (b) for 60, 90, 
and 480 ppm samples of HMW polymer. 

Compelling evidence that hydrodynamic interactions are important also comes from brow- 
nian dynamic simulations of A-DNA, a biopolymer similar in size and nature to our HMW 
sample [8], in shear and elongational flows [29] which demonstrate that for — above 0.1 in- 
cluding hydrodynamic interactions significantly affects calculated properties. The excellent 
agreement of the ~ 0.1c* threshold found in the simulation of (bio)polymers where hydro- 
dynamic interactions become important and the peaking of fluctuations here at « 0.1c* 
supports the hypothesis that hydrodynamic interactions are limiting growth in fluctuations. 
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We note that below the predicted threshold, above which hydrodynamic interactions 
have significant affect on bulk properties, our measure has approximate (square) power-law 
scaling as expected for the variance given linear growth in amplitude, which is in turn is 
expected if polymer coils are independent. The approximate power-scaling occurs with a 
concave rounding versus the strictly linear curve predicted on a log-log plot. This round- 
ing (breaking of power-law scaling) suggests subtle polymer-polymer interactions even in 
quite dilute solutions; recent systematic and careful rheological studies have demonstrated 
[HTj polymer-polymer interactions can exist for concentrations as low as 1% of the overlap 
concentration in elogational flows. 

The lack of true power-law scaling at low concentrations, the mild decrease in fluctua- 
tions above a threshold, and the correspondence of this observed threshold with the predicted 
threshold for hydrodynamic interactions in (hard) colloids and (soft) polymers indicate hy- 
drodynamic interactions are important in nominally dilute solutions and limit fluctuation 
growth. 
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FIG. 6: Fluctuations at various concentrations. The fluctuations in normalized PMT voltage for 
HMW polymer at different concentrations and 0.2 kV. It can be seen that fluctuations increase 
with concentration but rapidly peak and plateau well below the overlap concentration (« 0.1c*; 
c* ~ 300 here) . Modeling the polymeric solution as a hard colloidal suspension predicts a transition 
in behavior at c ~ 40-60 ppm (grey band) due to molecular crowding. Simulations of a similar 
polymer [29J find a threshold of 0.1c* (dotted line) above which hydrodynamic interactions between 
polymer coils significantly affects properties. The good agreement of both the colloidal model and 
the simulations with the observed peaking of fluctuation growth suggests hydrodynamic interactions 
and molecular crowding is the mechanism inhibiting fluctuations. 
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The relaxation time for our polymers is estimated with 

A Zimm « (1) 

Here Azimm is the Zimm relaxation time estimate [52] and r/ s the solvent viscosity that 
immerses the individual polymer coils of radius of gyration R g . The solvent frictionally slows 
down relaxation which is itself driven at a given thermal energy k^T which sets the level of 
agitation causing polymers to collapse into a random globular coil. 

It should be noted that in studies of viscoelastic fluids and their instabilities high viscosity 
solvents are most often used [T] , as this both increases the relaxation time due to the scaling 
with 77,3 and reduces the effects of diffusion which will tend to reduce contrast if dye is used 
for visualization. For EOF driven fluids velocity will scale inversely with rj s , the effect of 
increasing viscosity will therefore both increase relaxation time and reduce the achievable 
flow rates in such a manner that the effects will cancel leaving the De number unchanged. We 
initially attempted flows of high viscosity solvents (~ 65 % sugar and equivalent glycerine), 
however we found it difficult to work with such flows due to long start up times and low 
achievable velocities making it difficult to study the flows. Joule heating effects further 
complicate highly viscous flows. As low viscosity solvents are most often used for microfluidic 
applications, De is not expected to change in a significant manner with changing viscosity for 
EOF flows, and experimental studies were complicated by the use of high viscosity solvents 
we moved to our low viscosity solvent. 

Dynamic Light Scattering experiments (Brookhaven Instruments 200SM) for 120 ppm of 
the HMW sample gives a hydrodynamic radius of 296 nm while the LMW sample at 120 
ppm has a radius of 137 nm. Using the experimental relationship [31] between the radius of 
gyration R g and the hydrodynamic radius for polyacrylamide (i?h = 0.68i? g ; close to the 
Kirkwood-Riesman prediction |35j of R^ = 0.665i? g ) allows us to determine the relaxation 
times of 0.025 and 0.0025 s for the HMW and LMW samples, respectively. Using the Mark- 
Houwink parameters [35] for PAAm in pure water gives relaxation times of 0.028 and 0.0029 
s. The similarity in predicted relaxation times indicates that the addition of 20% MeOH 
has only a mild effect on the radius of gyration. We take the relaxation times as 0.03 and 
0.003s. 
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The viscoelastic instability condition (see section 2.8 in [T], see also [37]) is given by 

Xv XzimmV Xzimm^EOFE 1 

De = ~ = > -. 

-^Char -^Char -^Char 2 

That is, electric fields give rise to extensional flows [3], with significant coil stretching, 
when De > 1/2 pQ |37J for dilute solutions; note that for melts a threshold of 0.5 is also 
reported, which indicates that this approximate threshold value may be universal [3E]- For 
the HMW sample we calculate V cr ~ 0.015 kV for our system [53J. Above this threshold 
instability is observed, however below w 0.05 kV bias it becomes difficult to obtain data, 
due to diffusion rates becoming non-negligible for the low flow speeds, and the fluctuation 
in signal approaching the noise floor. 

To better observe the transition to instability the LMW sample is used at a concentra- 
tion of 120 ppm. The predicted critical driving voltage for the LMW sample is over an 
order of magnitude larger, V CI ~ 0.28 kV (corresponding to a critical velocity of roughly 2 
mm/s), than for the HMW sample. The threshold observed in Fig. [7J which occurs near 
the predicted value, supports the hypothesis that polymer stretching drives the instability 
[54] . The discrepancy between the observed and predicted threshold we attribute to both 
the relaxation time estimate and to existence of molecular weight distribution in polymer 
samples. Additionally, the theoretical critical Deborah number of 0.5 is approximate, and 
in direct experimental observation of DNA polymers in extensional flows [41 J at a Debo- 
rah number of [55] ~ 0.4 significant stretching occurred (the latter result would reduce the 
estimated V cr here to ~ 0.23 kV, a velocity of 1.6 mm/s). Within the limited precision of 
predictions our experimentally observed instability threshold agrees well with theory. 

Two alternative instability mechanisms must be discussed. Conductivity gradients be- 
tween co-flowing streams can give rise to electrohydrodynamic (EHD) instabilities when 
driven by electric fields [42]. The possibility of EHD instability is ruled out, as the (oth- 
erwise identical) polymer free solutions do not show instability. It is also conceivable that 
dynamic polymer-wall adhesion [19] of polymer coils coating the walls could lead to a time 
varying /^eof, driving instabilities. We reject this as underlying the instability, as PAAm 
does not adhere to glass [T5]. This rejection is also supported by the similarity between 
/zeof measurements of polymer-free and dilute (120 ppm) HMW polymer solutions which 
indicates no significant modification at the wall occurs. 
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FIG. 7: Fluctuation growth transition. A 120 ppm LMW sample (circles) with a predicted elastic 
instability at v CT ~ 2 mm/s experimentally demonstrates an instability onset at ~ 1.5 mm/s. Above 
~ 1.5 mm/s fluctuations grow dramatically, as compared to a gentle increase below the observed 
threshold; polymer free solutions (squares) display no significant fluctuations. 

IV. FORMATION OF GLOBULES 

We found that at large applied voltages globules would form which appeared to be gelled 
rafts of polymer (see Fig. [8]). The observed structures retained their form while the globule 
remained in the field of view. These globules only formed in the HMW polymer sample, 
at higher voltages (> 0.5 kV) and c > 60 ppm; at lower voltages, or for the LMW sample 
at 120 ppm, no indication of gel formation was observed. At the voltages reported in the 
present study (< 0.4 kV) no visual evidence of gel formation was observed in the HMW 
sample, and at no voltage was gellation observed for the LMW sample. 

The mechanism underlying gel formation is likely thermally activated cross-linking at the 
(Pt) electrode; electrogeneration of polymers at metal surfaces is an established synthesis 
technique which requires high current densities [33! • The possibility of contact glow elec- 
trolysis is discounted as current-voltage characteristics are linear both in the presence and 
absence of globules, which is not consistent with contact glow electrolysis where linearity 
is only observed pre-ignition of the plasma [33] • Dip coating [35] the electrodes and inten- 
tionally hard drying polymer films in place could either prevent flow, or resulted in globule 
formation at lower voltages (0.2 kV), and removal of films prevented globule formation at 
lower voltages. This result supports a high current density/thermal route to globule forma- 
tion where pin hole or other imperfections in hard dried films would lead to a high local 
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current density and subsequent heating effects. It is also known that polyacrylamide will 
spontaneously adhere to metal and form a partial film this partial film will result in 
higher current densities through the non-coated Pt surface relative to an pristine surface 
and a mildly higher surface coverage appears to arise for higher molecular weight samples 
[35] . This affinity for coating electrodes requires careful cleaning between runs [56], while a 
dependence of coverage on molecular weight could account for no observed globules arising 
in the LMW samples for voltages observed to give rise to globules in the HMW samples. 
Spontaneous coating of electrodes is a rapid effect [IB], on the order of 10's of seconds, and 
experimental setup times ensure that the coating process is complete before experiments 
commence. 



FIG. 8: Globule formation. Globules are formed at larger applied voltages (>0.5 kV). Seen here 
are globules formed at 1 kV potential for a 120 ppm HMW sample. It was found that at larger 
applied voltages globules would form in the fluid with dimensions on the order of the channel, 
presumably due to cross linking of PA Am at the electrodes. While this effect may be of interest 
in making microgels |47j on demand in microfluidics it is not explored further here. 



V. CONCLUSION 

In summary, extensional instabilities are excited in polymer solutions electro-osmoticaly 
pumped through a 2:1 microchannel constriction. Polymer-free solutions display creeping 
laminar flow, while addition of high molecular weight polymer leads to instabilities above 
a critical flow rate corresponding to the viscoelastic instability condition for extensional 
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flows. Electo-osmotic pumping is increasingly relevant as dimensions scale down and vis- 
coelastic effects become more pronounced for small dimensions due to the inverse scaling 
of the Deborah number with channel size (see equation |2|. We observe a dependence in 
measured fluctuations on polymer concentration which rapidly plateaus well below the over- 
lap concentration, occurring near the threshold where hydrodynamic interactions between 
polymer coils become important. At higher voltages microgels were formed at the electrodes 
and moved downstream. Our experimental results demonstrate instabilities can be readily 
excited, even in microscaled devices where laminar flow normally predominates. Electro- 
osmotic flow can be made unstable intentionally through the addition of small amounts of 
chemically inert polymer allowing the possibility for enhanced mixing. Our results also sug- 
gest caution when electro-osmically pumping biological [48, 49 j , or otherwise viscoelastic, 
samples where instabilities may arise under conditions typically assumed to be laminar. 
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